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Kinematic structure in the young a Orionis association 
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ABSTRACT 

We have used precise radial velocity measurements for a large number of candidate 
low-mass stars and brown dwarfs, to show that the young a Ori "cluster" consists of 
two spatially superimposed components which are kinematically separated by 7kms~ 1 
in radial velocity, and which have different mean ages. We examine the relationship 
of these two kinematic groups to other populations in the Orion OBI association and 
briefly discuss the consequence of mixed age samples for ongoing investigations of the 
formation and evolution of low-mass objects in this much-observed region. 

Key words: stars: low-mass, brown dwarfs - stars: pre-main-sequence - open clusters 
and associations: individual: a Orionis. 



1 INTRODUCTION 

Since their discovery (Walter et al. 1997), the low-mass ob- 
jects surrounding the young star a Ori (a quintuple system 
with an 09.5V primary and four spatially resolved, early- 
type components) have become one of the most intensively 
studied populations of young stars and brown dwarfs. They 
are reasonably close, so the young, intrinsically luminous, 
low-mass objects are bright. Interpretation of results is fur- 
ther aided by a high galactic latitude and low extinction 
(Brown, de Geus & de Zeeuw 1994). 

Examples of recent investigations include the search for 
very low-mass stars, brown dwarfs and even "free floating 
planets" and determination of the initial mass function at 
low masses (e.g. Zapatero Osorio et al. 2000, 2002a; Bejar 
et al. 2001; Sherry, Walter & Wolk 2004); using the asso- 
ciation to constrain the timescales for disc dissipation and 
the cessation of accretion in low-mass objects (e.g. Barrado 
y Navascues et al. 2003; Oliveira et al. 2004); searching for 
analogues of classical T-Tauri stars among young substellar 
objects (e.g. Jayawardhana et al. 2003); studying the evo- 
lution of magnetic activity and angular momentum in low- 
mass stars and brown dwarfs (e.g. Franciosini, Pallavicini & 
Sanz-Forcada 2006); and, investigating variability and possi- 
ble age spreads in young star forming regions (e.g. Caballero 
et al. 2004; Burningham et al. 2005a). 

Most of these investigations rely on the assumption of 
either a distance or age for the low-mass objects in the as- 
sociation. Furthermore, an implicit assumption often made, 
is that the studied objects are coeval or at least have an age 
spread that is modest compared with the cluster age. Previ- 
ous papers have either adopted the Hipparcos distance to a 
Ori (352^gg 6 pc - Perryman et al. 1997) and deduced an age 



of ~ 5 Myr from low-mass model isochrones (e.g. Oliveira 
et al. 2004), or they have assumed that a Ori is part of the 
larger Orion OBlb association that has a mean Hipparcos 
distance of ~ 440 pc (Brown, Walter & Blauuw 1998) and 
hence deduced an age of ~ 2.5 Myr (e.g. Sherry et al. 2004). 

In this letter we present the initial results of a precise 
radial velocity survey of low-mass stars and brown dwarfs 
around a Ori. We show that the "cluster" of low mass ob- 
ject around a Ori actually consists of two, kinematically 
distinct, spatially superimposed populations, with different 
mean ages. 



2 OBSERVATIONS AND ANALYSIS 

Multi-object spectroscopy was performed at three epochs 
for four separate fields (see Fig. around a Ori during 
the period 10 November 2005 to 13 January 2006. A total of 
148 targets were observed at each epoch using the FLAMES 
instrument mounted on the VLT-Kueyen (UT2) 8.2-m tele- 
scope. The Medusa fibre system fed the Giraffe spectrograph 
and a 316 lines mm -1 grating produced spectra with a re- 
solving power of 16000, covering the wavelength range 8080- 
8630A at 0.14A per pixel. Targets with 14.1 < I < 19.3 were 
selected from the I versus R — I colour-magnitude diagram 
(CMD) described by Kenyon et al. (2005), corresponding to 
masses of 0.5 > M/Mq > 0.03 (for an age of 3 Myr and 
distance of 440 pc - Baraffe et al. 1998, see Fig.gJ. 

A single exposure of 2750 s was obtained at each epoch, 
yielding spectra with signal-to-noise ratios (SNR) of 5-100 
per pixel. One or two bright F/G stars were simultaneously 
observed in each field using the FLAMES fibre feed to the 
UVES spectrograph. These gave very high SNR spectra at 
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Figure 1. A histogram of heliocentric radial velocities for the 
filtered subsample described in the text. Two kinematic groups 
with small internal velocity dispersions are identified. 



a resolving power of 47000, which covered the same wave- 
length range and allowed us to model telluric absorption. 

Spectra were reduced and extracted using the FLAMES 
UVES and Giraffe pipeline software. Giraffe wavelength cal- 
ibration was initially provided by thorium-argon lamp expo- 
sures and then improved by reference to telluric absorption 
lines in each extracted spectrum. Sky was subtracted by 
forming a clipped mean sky spectrum from the many fibers 
allocated to sky in each configuration. Telluric correction 
was provided by modelling telluric features in the simul- 
taneous UVES spectra using a 6-layer atmospheric model 
(Nicholls 1988) and the HITRAN molecular database (Roth- 
man et al. 2005). This synthetic telluric spectrum was con- 
volved with a gaussian profile and divided into the lower 
resolution Giraffe spectra. This removed all visible traces of 
telluric contamination. 

Heliocentric radial velocities (RVs) were found by cross- 
correlation with the template star HD 34055 (M6V) using 
the wavelength range 8061-8535A, but masking out regions 
around sky emission lines. Internal uncertainties in the RVs 
are < lkms -1 for most targets. The zero point was set 
by cross-correlating HD 34055 with similar UVES spectra of 
the M-dwarfs Gl 402, Gl 406 and Gl 876, which have precisely 
measured heliocentric RVs (see Bailer- Jones 2004). The he- 
liocentric RV of HD 34055 was +44.2kms _1 , with a range 
when using the three standards that leads to an external 
error estimate of ±0.5 km s -1 . 

Equivalent widths (EWs) for the gravity sensitive Nal 
(8183, 8195A) features were found by integrating in a 
±120kms range around the line centres (corrected for the 
RV). For the following analysis steps we then excluded: ob- 
jects with a significantly variable RV over the three epochs 
(the topic of binarity in the sample will be tackled in a sub- 
sequent paper); all objects with an EW[Na] < lA (indicative 
of a giant spectral class - see below); and all objects with 
a constant RV outside the range 10 < V r < 40kms _1 . A 
histogram of the weighted mean RVs for the remaining 117 
objects is shown in Fig. The RV distribution is clearly 
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Figure 2. The spatial distribution of the filtered subsample de- 
scribed in the text. The two kinematic groups are indicated by 
triangles (group 1) or circles (group 2). Small crosses are targets 
in neither kinematic group (wrong velocity, variable velocity or 
low EW[Na] - see text) and the four fields observed with the 
Giraffe spectrograph are indicated with large circles. 
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Figure 3. The total equivalent width of the Nal (8183, 8195A) 
doublet as a function of R — I colour. Members of the two kine- 
matic groups are indicated (triangles - group 1, circles - group 
2) and compared with data for field M dwarfs obtained with sim- 
ilar spectral resolution (Xu 1991; Montes & Martin et al. 1998; 
Kenyon et al. 2005). Small crosses mark targets in neither kine- 
matic group. 



bimodal and we refer to these two populations as kinematic 
groups 1 (20 < V r < 27kms _1 ) and 2 (27 < V r < ^^kms' 1 ) 
respectively. 

Figure [5] shows the spatial distribution of objects in the 
two kinematic groups. Representatives of both groups are 
found in all 4 Giraffe fields (and at all epochs) , but the ratio 
of group 2 to group 1 objects becomes smaller towards the 
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Figure 4. The colour magnitude diagram for the two kinematic 
groups (triangles — group 1, circles — group 2). Small crosses 
mark spectroscopic targets in neither kinematic group. Isochrones 
overplotted from Baraffe et al. (1998) are for ages/distances of 
3Myr/439pc and 10Myr/326pc. Reddening and extinction cor- 
responding to E(R — I) = 0.036 has been applied. 

north (e.g. this ratio is 10/15 in the north-west Giraffe field 
compared with 46/4 in the south-west field). Note that the 
photometric selection criteria for each of these fields were 
identical. The larger total number of targets in the south- 
west field simply reflects the larger number of photometric 
candidates near a Ori - a fact previously used to support 
the existence of the "cr Ori cluster" (e.g. Sherry et al. 2004) . 

Figure |3| shows how EW[Na] behaves with colour for 
each group. This feature strengthens considerably with grav- 
ity in cool stars (e.g. Schiavon et al. 1997). At young ages, 
low-mass stars and brown dwarfs contract almost isother- 
mally, resulting in a factor of ~ 10 increase in gravity be- 
tween 1 and 100 Myr at a similar colour (Baraffe et al. 1998) . 
Objects in kinematic groups 1 and 2 have a weaker EW[Na] 
than field M dwarfs of similar colour but stronger than giants 
where EW[Na] ^ lA is found. Hence they have gravities in- 
termediate between dwarfs and giants as expected for young 
low-mass objects. Figure |3] also shows that on average the 
objects in group 2 have lower gravities than those in group 1, 
although there is considerable overlap in the distributions. 
Formal EW uncertainties are less than 0.1 A for most of our 
targets and so cannot account for the overlap. Qualitatively, 
this indicates that group 1 is older (on average) than group 
2, but the lack of any detailed EW-age calibration prevent 
us from making a quantitative estimate of the age difference 
or any age spread within each group. 

Figure0]shows the I versus R—I CMD for the two kine- 
matic groups as well as for all the other spectroscopic targets 
to illustrate the target selection region. No obvious difference 
in distribution is seen. Of course all the targets were selected 
using the same photometric criteria, but Burningham et al. 
(2005b) have shown that few young stars exist outside of this 
range near a Ori. Therefore this diagram demonstrates that 
it is impossible to separate the two kinematic populations 
on the basis of their colours and magnitudes. 



3 DISCUSSION 

3.1 The link with Orion OBI 

The RV distribution shown in Fig. is clearly bimodal. 
Kinematic groups 1 and 2 have weighted mean RVs of 
(23.8±0.2) kms" 1 and (31.0±0.1) kms" 1 respectively (with 
an external error of ±0.5 kms -1 ). The standard deviations 
about these mean velocities are 1.1 kms -1 and 1.3 kms -1 , 
suggesting coherent kinematic groups with internal disper- 
sions similar to those found in young open clusters. Previ- 
ous RV investigations of stars around a Ori had insufficient 
precision (a ~ 5 kms -1 ) to properly resolve these two com- 
ponents (e.g. Walter et al. 1997; Kenyon et al. 2005). From 
Figs. EJand|3]we deduce that: (i) the two kinematic com- 
ponents have a different spatial distribution; (ii) on average, 
group 1 is older than group 2 but by an uncertain amount 
and they may overlap; and (iii) the two groups occupy sim- 
ilar positions in the CMD. 

The early- type stars of the Orion OBI association were 
originally split into four subgroups by Blaauw (1964), pri- 
marily on the grounds of position. The two subgroups rel- 
evant to this study are the la and lb subgroups - OBla 
is an older dispersed population mainly to the north and 
north-west of the belt stars, while OBlb consists of the belt 
region including the belt stars and a Ori. Brown et al. (1994, 
1998) have reviewed the subgroup properties using new pho- 
tometry and Hipparcos parallaxes for the O and B stars. 
They find the subgroups can also be separated on the basis 
of age and distance, obtaining ages of 11.4 ± 1.9 Myr and 
1.7 ± 1.1 Myr and distances of 326 ± 16 pc and 439 ± 33 pc 
for OBla and OBlb respectively. The two subgroups can- 
not be distinguished kinematically. Proper motions are very 
small in Orion OBI because the association is moving al- 
most radially away from the Sun. Average radial velocities 
of (23.8 ± 0.7) kms -1 and (23.1 ± 1.4) kms -1 are quoted 
for subgroups OBla and OBlb by Morrell & Levato (1991) 
based on measurements of large numbers of early-type stars. 
Bricefio et al. (2005) have searched for low-mass PMS stars 
in the Orion OBla and OBlb regions using variability as an 
identifying characteristic and spectroscopy for confirmation. 
Using Blaauw's spatial criteria to define the subgroups and 
assuming the distances quoted above, they compare their 
data with low-mass isochrones, finding ages of 7-10 Myr and 
4-6 Myr for OBla and OBlb respectively. The difference in 
age is supported by the frequency with which accretion discs 
are found in the two populations. 

All of our targets are spatially within Blaauw's OBlb 
region, but given that subgroup OBla is in the foreground 
and widely dispersed, there seems little reason to suppose 
there are not significant numbers of older OBla members in 
the line of sight towards OBlb (e.g. see Fig. 9 in Sherry et 
al. 2004) . Hence, based on the similarity of their mean RVs, 
kinematic group 1 could be identified with either the OBla 
or OBlb subgroups. Conversely, kinematic group 2 cannot 
be identified with any of the Orion OBI subgroups and has a 
mean RV which is similar to the RV of cr Ori itself, a Ori AB 
has V r — 29.9 ± 1.6km s -1 according to Barbier-Brossat & 
Fignon (2000), who analyse all pre-1990 RV data. Morrell & 
Levato (1991) obtain 27±4kms -1 for a Ori E, but caution 
that it is probably variable. Given this kinematic evidence 
it is hard to see (other than by a simple positional criterion) 
how cr Ori could be securely classed as part of Orion OBlb 
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and therefore it should not be assumed that either a Ori or 
the group 2 objects around it have a distance and age in 
common with Orion OBlb. 

Figure 0] shows two isochrones taken from the models 
of Baraffe et al. (1998) and converted into the observational 
plane using the empirical colour-effective temperature rela- 
tionship and bolometric corrections described in Jeffries et 
al. (2004). These isochrones are approriate for our current 
knowledge of the age and distance of the Orion OB la and 
OBlb subgroups. Our targets, along with Orion OBla and 
OBlb, lie close to a locus (or more correctly a band) in 
the age-distance plane that results in young coeval popula- 
tions appearing in roughly the same position in the colour- 
magnitude (and Hertzsprung-Russell) diagram. 

The RV evidence, the strong spatial concentration of 
low-mass PMS stars and the low ratio of group 1 to group 
2 objects immediately around a Ori strongly support the 
identification of group 2 as the "cr Ori cluster". The Hip- 
parcos parallax of a Ori only constrains the distance to this 
cluster to be 270-520 pc, with a corresponding age from low- 
mass isochrones of between 15-2 Myr. However, if kinematic 
group 1 is older on average than group 2 (as suggested by 
Fig. [31 and, identified with Orion OBla at ~ 10 Myr, then 
this would favour an age of a few Myr for group 2 and a dis- 
tance similar to that of Orion OBlb (~ 440 pc). This would 
be consistent with a likely main-sequence lifetime of < 7 Myr 
for a Ori A (09.5V, T eS ~ 3.0 x 10 4 K - Schaller et al. 
1992). Alternatively, if group 1 is dominated by members of 
the more youthful Orion OBlb, then a Ori and the group 2 
objects that surround it would have to be even younger and 
further away. 

There are two ways to further constrain this prob- 
lem. One would be to improve the parallax measurement 
to a Ori and the second would be to investigate the two 
kinematic groups with distance-independent age indica- 
tors. Of the latter there are two possibilities. Photospheric 
lithium is rapidly depleted in low-mass stars beginning at 
an age which is mass (and hence temperature dependent). 
Zapatero-Osorio et al. (2002b) found that all their a Ori as- 
sociation targets across a wide mass range had Li consistent 
with no depletion, leading to an upper age limit of 8 Myr. 
Unfortunately their targets do not have precise RVs and 
few are contained within our survey, so we do not know to 
which kinematic group they belong. Instead we have cross- 
correlated members of kinematic groups 1 and 2 with the 
targets of Kenyon et al. (2005), finding 39 correlations (10 
in group 1 and 29 in group 2). Li-rich stars are found at 
R — I 1.6 in both groups (warmer stars were not observed 
by Kenyon et al.), corresponding to an effective temperature 
of ~ 3400 if. Using the models of Baraffe et al. (1998), this 
places an upper limit of < 15 Myr to the age of both groups. 
A similar limit is found using other evolutionary models. De- 
tailed Li spectroscopy of warmer members of groups 1 and 
2 might be able to improve these constraints. 

A second empirical method is to compare the fraction of 
objects with circumstellar material with those in other clus- 
ters and assume that disc dissipation occurs on a universal 
timescale. We have correlated objects in kinematic groups 1 
and 2 with objects observed at JHK S L' by Oliveira et al. 
(2004, 2006). We find 29 correlations and that 1/4 members 
of group 1 and 10/25 members of group 2 have a significant 
K s — L' excess (as defined in those papers). There are in- 



sufficient statistics to comment on the relative ages of the 
two groups, but a 40 per cent disc frequency for group 2 is 
consistent with an age of 3-4 Myr according to Haisch, Lada 
& Lada (2001). 

In summary, our favoured, although not necessarily 
unique, interpretation is that kinematic group 1 should be 
identified with Orion OBla at an age of ~ 10 Myr and at 
a distance of ~ 330 pc. The spread about a single isochrone 
could be due to the presence of some younger Orion OBlb 
members and perhaps a large spread of distances within the 
widely dispersed OBla association. Group 2 would then have 
an age of about 3 Myr, be clustered around a Ori and be at 
a similar distance to Orion OBlb (~ 440 pc). However, the 
"cr Ori cluster" appears to be kinematically distinct from 
Orion OBlb. 

3.2 Consequences 

Studies of low-mass stars and brown dwarfs in the a Ori 
region must account for the possibility of a population with 
mixed age and distance. These populations cannot be distin- 
guished in the CMD and the only way to reliably attribute 
stars to one group or the other is via RV measurements 
with a precision better than a few kms -1 . The populations 
do show some spatial segregation. It appears that samples 
restricted to within 10 arcminutes of a Ori would suffer con- 
tamination from the older population at the level of only 
~ 10 percent. However, in the quest for larger and more com- 
plete samples, many studies in the literature have travelled 
beyond this boundary, into regions (e.g. > 15 arcminutes to 
the north of a Ori) with roughly equal numbers from both 
populations. How much of a problem would older interlop- 
ers be? There is insufficient space in this letter to explore in 
detail the effects of our discovery on all the previous studies 
of the a Ori region but we make some general points below. 

The effects on mass determination depend on which 
technique is used. A simple mass- magnitude relationship ob- 
tained from models (e.g. Baraffe et al. 1998) is not greatly af- 
fected by assignment of a star to one kinematic group or the 
other - the differences in age and distance almost compen- 
sate. Masses for an age of 3 Myr and distance of 440 pc are 
larger than those deduced for an age of 10 Myr and distance 
of 330 pc by about 6 per cent at 0.5 Mq, by 30 per cent at 
0.06 A'Iq and are almost the same at 0.02 Mq. However, big- 
ger discrepancies would arise if masses were estimated from 
colours or spectral types because the evolutionary tracks are 
not vertical when using these parameters as abcissae in the 
Hertzsprung-Russell diagram. At a given colour (in the I 
versus I — J CMD for example) one might deduce that a 
star had 0.5 Mq by assuming an age of 3 Myr, but it would 
actually be 0.25 Mq if the age were 10 Myr (Baraffe et al. 
1998). This discrepancy reduces towards lower masses as the 
evolutionary tracks become more vertical in the CMD, be- 
coming a factor of 1.2 at O.O7M and ~ 1 below 0.04 Mq. 
There are also significant structural changes for low-mass 
stars on these timescales which would be important if using 
a <r Ori sample to study the evolution of angular momen- 
tum and rotation rates. Under the assumptions of negligible 
angular momentum loss between 3 and 10 Myr and homolo- 
gous contraction, objects thought to be at 3 Myr but which 
were actually at an age of 10 Myr, would have spun up by 
a factor of 2.2 at 0.5 Mq and 2.9 at 0.07 Mq . The inter- 
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lopers could therefore have a significantly different rotation 
distribution, perhaps appearing as a tail of short period, fast 
rotators. 

If the disc lifetime distribution is independent of envi- 
ronment then interlopers from an older population would 
have significantly fewer (accretion) discs and result in a 
downwardly biased estimate of the disc frequency in the 
a Ori cluster. As discussed in the previous section there 
are insufficient data to test this in detail yet. The disc fre- 
quency we determined for members of kinematic group 2 is 
higher than, but still consistent with, the 33 ± 6 per cent 
frequency determined for a larger sample of objects around 
a Ori which contains at least some objects from group 1. 

Finally, studies of the age spread and spatial distri- 
bution within the cluster must be aware of, and try to 
eliminate, the effects caused by including two populations 
with probable differences in age and distance. For exam- 
ple, the radial profile of the cluster deduced by Sherry et 
al. (2004) includes photometrically selected cluster members 
more than 30 arcminutes to the north of a Ori, which are 
quite likely to belong to kinematic group 1. Excluding these 
would reduce the number of stars and brown dwarfs in the 
a Ori cluster and decrease the deduced cluster core radius. 



4 CONCLUSIONS 

The cluster of young objects situated towards a Ori con- 
sists of two kinematic components separated by ~ 7kms~ 
in RV. The first component (and larger in our sample) is 
concentrated around a Ori and shares a common RV with 
it. The second component increases in concentration towards 
the north and has a RV similar to the Orion OBla and OBlb 
sub-associations. On the basis of gravity-sensitive spectral 
features, the second component appears to be older on av- 
erage than the first, although considerable overlap is possi- 
ble. The two components cannot be separated in a colour- 
magnitude diagram, a degeneracy that implies that the older 
second component is closer. Our favoured interpretation is 
that the first component is the "er Ori cluster" , but that it 
should not necessarily be identified (in distance or age) with 
any part of the Orion OBI association. The second com- 
ponent consists of low-mass objects from either the Orion 
OBla or OBlb sub-associations or more probably, a mix- 
ture of the two. The presence of multiple populations with 
different ages complicates the study of low-mass objects in 
this region. 
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